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Abstract. Galactic nuclei are unique laboratories for the study of processes connected 
with the accretion of gas onto supermassive black holes. At the same time, they represent 
challenging environments from the point of view of stellar dynamics due to their extreme 
densities and masses involved. There is a growing evidence about the importance of the 
mutual interaction of stars with gas in galactic nuclei. Gas rich environment may lead to 
stellar formation which, on the other hand, may regulate accretion onto the central mass. 
Gas in the form of massive torus or accretion disc further influences stellar dynamics in the 
central parsec either via gravitational or hydrodynamical interaction. Eccentricity oscilla- 
tions on one hand and energy dissipation on the other hand lead to increased rate of infall 
of stars into the supermassive black hole. Last, but not least, processes related to the stellar 
dynamics may be detectable with forthcoming gravitational waves detectors. 
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1. Introduction 

Standard model of active galactic nuclei con- 
sists of a supermassive black hole (SMBH) 
of mass Mbh surrounded by an accretion disc 
which extends to ~ 10 4 - 10 5 R g ; R g = 
GM Bli /c 2 . Further away, at r > 10 6 R S , there is 
assumed to be an obscuring massive molecular 
torus. Within the radius of ~ 10 6 - 10 7 /? g , there 
is a dense stellar cusp, of the total mass com- 
parable to Mbh- In spite of its extreme density 
which may exceed 10 8 stars per cubic parsec 
in the centre, the two-body relaxation time in 
this region is of the order of Gyrs. Hence, on 
the orbital (crossing) time scale, the stars fol- 
low nearly Keplerian orbits in the gravitational 
field of the SMBH. Interaction of stars with the 



axisymmetric gaseous structures leads to a sec- 
ular orbit evolution. Its characteristic time may 
be well below the relaxation time and, there- 
fore, it can lead to observable effects. 

In this contribution we briefly introduce 
various modes of the mutual interaction of the 
stars and gas in galactic nuclei and discuss pos- 
sible consequences. In Section [4] we present 
our results within the context of the nucleus of 
the Milky Way. 

2. Gravitational interaction 

Gravitational field of the SMBH as well as the 
mean field of the star cluster are assumed to be 
spherically symmetric, which implies conser- 
vation of the vector of the angular momentum. 
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This is no longer true once a non-spherical 
perturbation to the gravitational field is intro- 
duced. We will concentrate on the case of an 
additional component of the gravitational field 
due to the accretion disc or molecular torus 
for which we assume an axial symmetry. Then, 
one component of the angular momentum vec- 
tor, L z , parallel to the symmetry axis will be 
conserved, but the eccentricity may change. 

This process, sometimes referred to as 
Kozai oscillations, is well described by means 
of the perturbation Hamiltonian theory for an 
analogical (reduced) hiera rchical triple system 
dKozailll96d lLidovlll962l) . For the simplest 
case of perturbation potential due to an in- 
finitesimally narrow ring of mass and ra- 
dius R^, the equations of motion in the first ap- 
proximation read: 
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7k rj = ~ cos i [l + 4e 2 - 5e 2 cos 2 <j] , (4) 

where a, e, i, a> and Q are the mean orbital 
semi-major axis, eccentricity, inclination with 
respect to the plane of symmetry of the perturb- 
ing potential, argument of the pericentre and 
the longitude of the ascending node, respec- 
tively. We denote rj = Vl - e 2 and 
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which is the characteristic time of the evolu- 
tion of the orbital elements. Equations (Q]) - 0) 
have two integrals of motion, 

Ci=77sin/ and C2 = (5 sin 2 i sin 2 a> - 2)e 2 , (6) 

which implies that three orbital elements, e, i 
and to, change periodically with equal period 
and phase. For low values of C\, the ampli- 
tude of the oscillations of eccentricity may, 
in terms of rj, reach several orders of magni- 
tude. Therefore, we expect enhanced rate of 
events connected with extreme eccentricities, 
e.g. tidal disruptions of stars. 
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Fig. 1. Size of the loss cone of tidal disruptions as a 
fraction of the phase space. Lower dotted line repre- 
sents size of the loss cone for an unperturbed case. 
Upper dotted line corresponds to an estimate for the 
case of perturbation due to the ring. Model repre- 
sented by the solid line contains the 1 .PN correction 
in addition. Dashed lines stand for complex models, 
including the gravity of the stellar cusp. 



In the case of the galactic nuclei, there are, 
however, two additional perturbations to the 
dynamics of individual stars, that tend to di- 
minish the amplitude of the eccentricity os- 
cillations. It is the mean potential of the stel- 
lar cusp and the relativistic pericentre advance. 
For some sets of parameters, both of them 
lead to substantial increase of |dw/df|, i.e. to 
the shortening of the period of the oscilla- 
tions and also to damping of their am plitude 
dBlaes et al. Il2002tllvanov et al. 1 12005b. 



In Karas&Subrl ( 2007 ) we have com- 



pared the volume of the loss cone of the cen- 
tral star cluster occupied by stars that reach 
the SMBH with in the tidal radius, R, . We 
have considered a lBahcall& Woif1 dl976) pro- 
file of the cusp and various disturbing com- 
ponents to the central Keplerian potential, i.e. 
the axially symmetric ring or disc, first post- 
Newtonian correction to the central potential 
and the spherically symmetric mean potential 
of the stellar cusp. The results are presented in 
Fig. [U When considering only the ring as a per- 
turbation, the loss cone is more than two orders 
of magnitude larger in comparison to the un- 
perturbed case. Relativistic pericentre advance 
damps the Kozai oscillations for a < a t , 



J ~ 32 p6 r>2 n-l ,,-2 

a t « -5- K d K g K t fx . 
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where we denote ji = M^/Mbh- This trans- 
forms into the break of the solid line in Fig.Q]at 
Mbh ~ 1O 6 M . The pericentre shift caused by 
the stellar cusp, mass of which was set equal to 
Mbh in Fig. Q] leads to a systematic decrease 
of the loss cone over the whole range of the 
black hole mass. We may conclude that ax- 
isymmetric perturbation to the potential of the 
SMBH may increase the number of the tidally 
disrupted stars by a factor > 10. This number 
may by further increased if the source of the 
axisymmetric component of th e potential ex- 
tends over a large range of radii dKaras & Subr I 
120071) . 

3. Hydrodynamical interaction 

The dynamical processes described in the pre- 
vious Section do not (in the first approxima- 
tion) change energy of individual stars. On 
the other hand, repetitive passages of the stars 

through the accretion disc lead to dissipa- 

i it — L — I 

tion of their kinetic energy (Sveret al. 1991; 
Vokrouhlickv&KarasI 119981; ISubr&KarasI 
1999). Being highly supersonic, the dissipa- 
tional efficiency of the star-disc collisions is 
assumed to be proportional to the physical 
cross-section of the star and the column density 
of the accretion flow. For a standard g a s pres- 
sure dominated Shak ura & Sunvaevl (119731) 
accretion disc the characteristic ti me of the or- 
bital decay can be estimated as ( Sub r et aTl 
120041) : 

f -^ M8 (|)fe) 9/4yr ' (8) 

where £„ = M*/Rl is the mean column den- 
sity of the star and M% = M B h/10 8 M q . Beside 
the shrinking of the semi-major axis this in- 
teraction leads to circularisation of the or- 
bit and its decay in to the plane of the disc 
(lKaras& Subr 11200 ll) . We distinguish two dif- 
ferent modes of subsequent migration towards 
the centre once the star gets embedded in 
the disc. If the star's Roche radius is larger 
than the half-thickness of the disc, it opens 
a permanent gap in it and follows its slow 
radial inflow. If the gap is not opened, the 
star excites density waves in the gas which 



leads to larger transfer of angular momen- 
tum and, consequently, fa s ter inward m igration 
(Lin&Papaloizou 1986; Ward 1986) on the 
time-scale of 

The inner region of the star cluster inter- 
acting with the accretion disc is assumed to be 
continuously supplied with new stars on the re- 
laxation time-scale: 

10! "''*»" ! (ilF < 10 » 

where ri(, = «o/(10 6 pc~ 2 ) is the stellar density. 
Condition t co \\ < f re i ax defines radius R out of the 
region which will be emptied due to the drag 
imposed by the disc upon the stars. The in- 
flow rate of stars towards t he centre through t he 
boundary of this region is dSubr et al. 112004 

fi-wrfe^**- <»> 

4. The Galactic Centre 

The closest galactic nucleus is that of the 
Milky Way. It harbours a SMBH of mass 
M BH ~ 3 5 x 10 6 M H dGenzel et al. 1 120031: 
iGhez et aT71 120031) which manifests itself as a 
radio source SgrA*. In spite of that current ac- 
tivity of SgrA*is highly sub-Eddington, there 
is an indirect evidence for an existence of 
massive accretion disc few million years ago: 
The near infrared observations of the central 
parsec have revealed a numerous population 
of young stars. Moreover, it has been found 
dLevin & Beloborodov 1 120031) that consider- 
able fraction of them forms a relatively thin 
disc rotating coherently around the SMBH. 
Fragmentation of a self-gravitating gaseous 
disc of mass > 1O 4 M is currently one of the 
most popular explanation of the origin of this 
young stellar structure. 

Stellar disc induces axisymmetric per- 
turbation to the potential of the SMBH. 
Therefore, we assume that it leads to the ec- 
centricity oscillations of late-type stars from 
the embedding spherical stellar cusp. In 
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iKaras&Subrl (120071) we have estimated that 
up to 100 stars (i.e. ~ 2% of the stars from the 
region with r < 0.3pc may have been tidally 
disrupted due to this process within the life- 
time of the disc. Another source of axisym- 
metric perturbation to the gravity of the central 
black hole is a circum-nuclear disc (CND), a 
molecular torus of mass > 0.1 Mbh and radius 
as 1 .6pc. In spite of its larger mass compared to 
the stellar disc, the CND is relatively less effec- 
tive in pushing the stars to extreme eccentrici- 
ties due to its larger radius. Smaller efficiency 
is, on the other hand, balanced by larger num- 
ber of stars that are influenced. Hence, another 
~ 100 stars from the region below 1.5pc may 
have been tidally disrupted due to the gravity 
of the CND. 

5. Conclusions 

Gas in the galactic nuclei can manifest itself 
not only by means of its own radiation. It forms 
complex environment together with the super- 
massive black hole and millions of stars. All 
these components mutually interact and influ- 
ence themselves. In this contribution we have 
focused on the effects of the gravitational and 
hydrodynamical influence of the gas in the 
form of accretion disc or torus upon the stars. 
We have discussed the effect of Kozai oscilla- 
tions and shown that they are likely to lead to 
enhanced rate of tidal stellar disruptions. Due 
to the relativistic pericentre advance, this pro- 
cess is strongly damped for Mbh <; 10 7 M o , 
while it is the most efficient for intermediate 
masses of the black holes, i.e. it can help them 
to grow from the IMBH to the SMBH stage. 

Passages of the stars through the gas dis- 
sipate their kinetic energy. This leads to their 
continuous inflow towards the black hole. The 
mass flow in the form of stars has to be small, 
compared to the accretion rate of the gas, oth- 
erwise, the gaseous structure would be de- 
stroyed. The individual stars may, however, be 
detectable in the final phase of their inspiral 
due to the emission of gravitational waves. 

Recent observations of the Galactic Centre 
indicate recent star formation at a distance 
of 10 5 - m b R g from the SMBH. This goes 
in line with previous theoretical considera- 



tions about fragmentation of the outer parts 
of the accretion dis cs due to self-gravity 
dCollin & Zahn 1119991) . In addition, recent nu- 
merical simulations of the fragmenting accre- 
tion flow suggest that star formation may in- 
hibit furth er inflow of the gas o nto the cen- 
tral mass (Nayakshi n et al. 1120071) . Hence, the 
feedback from stars should be considered in re- 
alistic models of AGNs. 
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